Nuclear matter from chiral low-momentum interactions 
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Nuclear matter calculations based on low-momentum interactions derived from chiral nucleon- 
nucleon and three-nucleon effective field theory interactions and fit only to few-body data predict 
realistic saturation properties with controlled uncertainties. This is promising for a unified descrip- 
tion of nuclei and to develop a universal density functional based on low-momentum interactions. 



PACS numbers: 21.65.-f, 21.30.-x, 21.60.Jz, 21.10.-Dr 

Nuclear forces saturate, so nuclei are self bound with 
roughly constant interior density. The Coulomb interac- 
tion drives heavier stable nuclei toward an imbalance of 
neutrons over protons and eventual instability, but the- 
orists can extrapolate a Coulomb-free, N = Z nucleus 
to arbitrary size. The uniform limit is called symmetric 
nuclear matter. For fifty years, an accurate prediction of 
nuclear matter starting from nuclear forces has been a 
theoretical milestone on the way to finite nuclei, but has 
proved to be an elusive target. Here we present the first 
nuclear matter calculations using soft Hamiltonians de- 
rived from chiral effective field theory interactions fit only 
to few-body (A ^ 4) data. We find realistic saturation 
properties with controlled uncertainties. 

Despite the long-term emphasis on the infinite uniform 
system, most advances in microscopic nuclear structure 
theory over the last decade have been through expanding 
the reach of few-body calculations. This has unambigu- 
ously established the quantitative role of three-nucleon 
forces (3NF) for light nuclei {A ^ 12) 0, 0. However, 
until now few-body fits have not sufficiently constrained 
3NF contributions at higher density such that nuclear 
matter calculations are predictive. Nuclear matter satu- 
ration is very delicate, with the binding energy resulting 
from cancellations of much larger potential and kinetic 
energy contributions. When a quantitative reproduc- 
tion of empirical saturation properties has been obtained, 
it was imposed by hand through adjusting short-range 
three-body forces (see, for example, Refs. [HHI). 

Progress for controlled nuclear matter calculations has 
long been hindered by the difficulty of the nuclear many- 
body problem when conventional nuclear potentials are 
used. The present calculations pull together several re- 
cent developments to overcome the hurdles: systematic 
starting Hamiltonians based on chiral effective field the- 
ory (EFT) 0, [||, renormalization group (RG) meth- 
ods @, H| to soften the short-range repulsion and short- 
range tensor components of the initial chiral interac- 
tions so that con verg ence of many-body calculations 
is vastly accelerated lfj, |llj, |l_2 1 , and an improved 3NF 



fitting procedure to the 4 He radius 0. In combination 
we obtain controlled theoretical uncertainties. 

Our results are summarized in Fig. [TJ which shows the 
energy per particle of symmetric matter as a function of 
Fermi momentum fcp, or the density p = 2kp/(3ir 2 ). A 
grey square representing the empirical saturation point 
is shown in each of the nuclear matter figures. Its 
boundaries reflect the ranges of nuclear matter satura- 
tion properties predicted by phenomenological Skyrme 
energy functionals that most accurately reproduce prop- 
erties of finite nuclei. Although this determination can- 
not be completely model independent, the value is gen- 
erally accepted for benchmarking infinite matter. In the 
future, calculations of the properties of finite nuclei will 
allow one to compare directly to experimental data. 

The calculations of Fig.[T]start from the N 3 LO nucleon- 
nucleon (NN) potential (EM 500 MeV) of Ref. [|. This 
NN potential is RG-evolved to low-momentum interac- 
tions Vi ow fc with a smooth 7i cxp = 4 regulator [§]. For 
each cutoff A, two couplings that determine the shorter- 
range parts of the N 2 LO 3NF [3] are fit to the 3 H 
binding energy and the 4 He matter radius using ex- 
act Faddeev and Faddeev-Yakubovsky methods as in 
Ref. [3. Our 3NF fit values are given in Table [0 We 
use the same 3NF regulator exp[— (p 2 + 3/4g 2 ) 2 /A| NF ] 
of Ref. [IH, but with a 3N cutoff A3NF that is al- 
lowed to vary independently of the NN cutoff. The 
shaded regions in Fig. [1] show the range of results for 
2.0 fm" 1 < A 3NF < 2.5 fnT 1 at fixed A = 2.0ml" 1 . 
Nuclear matter is calculated in three approximations: 
Hartree-Fock (left) and including approximate second- 
order (middle) and summing particle-particle-ladder con- 
tributions (right). The technical details are given in 
Ref. but we have improved the calculation to in- 

clude full momentum-dependent Hartree-Fock propaga- 
tors and (sub-leading) 3N double-exchange diagrams be- 
yond Hartree-Fock, however in a very approximate way. 
Further improvements are in progress [15j. 

The Hartree-Fock results show that nuclear matter is 
bound even at this simplest level. A calculation without 
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FIG. 1: (Color online) Nuclear matter energy per particle as a function of Fermi momentum &f at the Hartree-Fock level (left) 
and including second-order (middle) and particle-particle- ladder contributions (right), based on evolved N 3 LO NN potentials 
and 3NF fit to -Z?3 H and ?"4 He . Theoretical uncertainties are estimated by the NN (lines) and 3N (band) cutoff variations. 



approximations should be independent of the cutoffs, so 
the spread in Fig.[T]sets the scale for omitted many-body 
contributions. The second-order results show a dramatic 
narrowing of this spread, with predicted saturation con- 
sistent with the empirical range. The narrowing happens 
across the full density range. This is strong evidence that 
these encouraging results are not fortuitous, but that we 
have reached cutoff independence at the level of 1—2 MeV 
per particle. The 3NF fits to the 4 He radius improve the 
cutoff independence significantly compared to fitting to 
A = 3,4 energies only, see Fig. 6 in Ref. pj. For all 
cases, the compressibility K = 190 — 240 MeV (mainly 
from A3NF = 2.0 — 2.5 fm -1 ) is in the empirical range. 
To our knowledge, these are the first nuclear forces fit 
only to A ^ 4 nuclei that predict saturation reasonably. 

The particle-particle-ladder sum is little changed from 
second order except at the lowest densities shown. This 
is not surprising because at low density the presence of 
a two-body bound state necessitates a nonperturbative 
summation. We note that below saturation density, the 
ground state of nuclear matter is not a uniform system, 
but breaks into clusters (see, for example, Ref. [la]). 

In chiral EFT without explicit Deltas, 3N interactions 
start at N 2 LO p| and their contributions are given dia- 
grammatically by 




Ci,C 3 ,C4 CD Ce 



We assume that the Cj coefficients of the long-range two- 
pion-exchange part are not modified by the RG and take 
these values from Ref. [13]. At present, we rely on the 
N 2 LO 3NF as a truncated "basis" for low- momentum 3N 
interactions and determine the cd and ce couplings by 







SRG 


A or A/A 3N p [fm x ] 


CD Ce 


Cd 


ce 


1.8/2.0 


-0.0112 -0.2212 






2.0/2.0 


-0.3000 -0.2761 


— 1.023 


-0.3397 


2.0/2.5 


-2.000 -0.7564 


— 2.991 


-0.8797 


2.2/2.0 


-0.9000 -0.3673 






2.8/2.0 


-1.552 -0.4058 







TABLE I: Results for the cd and ce couplings fit to i?3 H = 
-8.482 MeV and r 4Ho = 1.95 - 1.96 fm for the NN/3N cutoffs 
used here. For Vi ow (SRG) interactions, the 3NF fits lead to 
Ei Hc = -28.22 . . . 28.45 MeV (-28.53 . . . 28.71 MeV). 



a fit to data for all cutoffs [3]. In the future, it will be 
possible to fully evolve three- and four-body forces start- 
ing from chiral EFT pi El- The fit values of Table Hare 
natural and the predicted 4 He binding energies are very 
reasonable. We have also verified that the resulting 3NF 
becomes perturbative in A — 3,4 (see also Ref.Jl4j) for 
the cutoffs used and generally for A3NF < A 15]. The 
sensitivity of many-body observables to uncertainties in 
the Ci coefficients was demonstrated recently in neutron 
matter calculations [20] . This raises the possibility of us- 
ing nuclear matter to constrain some of the Ci couplings. 

The evolution of the cutoff A to smaller values is ac- 
companied by a shift of physics. In particular, effects 
due to iterated tensor interactions, which peak in the 
relative momentum range k ~ 4fm _1 (and thus lead to 
saturation at too high density), are replaced by three- 
body contributions. The role of the 3NF for saturation is 
demonstrated in Fig. O The two pairs of curves show the 
difference between the nuclear matter results for NN-only 
and NN plus 3N interactions. It is evident that saturation 



3 




k p [fnf '] k p [fm' 1 ] 



FIG. 2: (Color online) Nuclear matter energy of Fig. [T]at the 
particle-particle-ladders level compared to NN-only results for 
two representative NN cutoffs and a fixed 3N cutoff. 



FIG. 3: (Color online) Nuclear matter energy of Fig. [T] at 
the particle-particle-ladders level comparing low-momentum 
Mow k with SRG-evolved chiral NN interactions. 



fm , which 



is driven by the 3NF [lfj. Even for A = 2 
is similar to the lower cutoffs in chiral EFT potentials, 
saturation is at too high density without the 3NF. Nev- 
ertheless, as in previous results [lfj, the 3N contributions 
and the cd-.ce fits are natural, and the same is expected 
for the ratio of three- to four-body contributions. 

The smooth RG evolution used in the results so far is 
not the only choice for low-momentum interactions. A 
recent development is the use of flow equations to evolve 
Hamiltonians, which we refer to as the Similarity Renor- 
malization Group (SRG) 2l|, l22j, 123J . The flow parameter 
A, which has dimensions of a momentum, is a measure of 
the degree of decoupling in momentum space. Few-body 
results for roughly the same value of SRG A and smooth 
Vjowfr. A have been remarkably similar (see, for example, 
Ref . ) . The analogous nuclear matter energies shown 
in Fig.[3]are also similar, which helps support the general 
nature of the 3NF fit. On the other hand, the difference 
of 2MeV per particle at saturation and above enlarges 
our theoretical uncertainty. 

The presented results are starting from a chiral EFT 
potential at a single order with one choice of EFT reg- 
ulator implementation and values. There are several al- 
ternatives available 0, @, [24]], which are largely phase- 
shift equivalent (the chi-square is not equally good up to 
-Eiab ~ 300 MeV). Universality for phase-shift equivalent 
chiral EFT potentials was shown for smooth-cutoff Viow fc 
interactions in Ref. [1] in the form of a collapse of the dif- 
ferent potentials to the same matrix elements. An anal- 
ogous collapse for N 3 LO potentials evolved by the SRG 
is shown in Fig.|4]for the ^^So channel, with a comparison 
to the corresponding Vj ow fc interactions. Similar results 
are found in the other channels. Based on this universal 
collapse, we do not anticipate large differences in nuclear 




EM 500 MeV 
EM 600 MeV 
EGM 550/600 MeV 
EGM 600/700 MeV I 
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FIG. 4: (Color online) Off-diagonal momentum-space ma- 
trix elements of different chiral N J LO potentials (EM 5] and 
EGM [(|) in the 1 So channel (upper panel) and after evolution 
to low-momentum Vi ow fe or SRG interactions (lower panel). 
A similar universal collapse is found for the low-momentum 
diagonal matrix elements and in other channels. 



matter when starting the evolution with different N 3 LO 
potentials, but a detailed comparison is forthcoming [FoT ] . 

The theoretical errors of our results arise from trunca- 
tions in the initial chiral EFT Hamiltonian, the approxi- 
mation of the 3NF, and the many-body approximations. 
We do not claim a chiral expansion for nuclear matter 
saturation, only that the hierarchy of potential energy 
contributions is maintained by the RG/SRG evolution. Is 
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the nuclear matter many-body calculation under control? 
Corrections to the present approximation include higher- 
order terms in the hole-line expansion and particle-hole 
corrections. While we have positive circumstantial ev- 
idence from cutoff independence that these corrections 
are small, we believe that an approach such as coupled 
cluster theory that can perform a high-level resummation 
is necessary for a robust validation. 

While nuclear matter has lost its status to light nu- 
clei as the first step to nuclear structure, it is still key 
as a step to heavier nuclei. Our results open the door 
to ab-initio density functional theory (DFT) based on 
expanding about nuclear matter (25| . This is analogous 
to the application of DFT in quantum chemistry and 
condensed matter starting with the uniform electron gas 
in local-density approximations and adding constrained 
derivative corrections. Phenomcnological energy func- 
tionals (such as Skyrme) for nuclei have impressive suc- 
cesses but lack a (quantitative) microscopic foundation 
based on nuclear forces and seem to have reached the 
limits of improvement with the current form of function- 
als [2^, [27]]. The theoretical errors of our results, while 
impressively small on the scale of the potential energy per 
particle, are far too large to be quantitatively competi- 
tive with existing functionals. However, there is the pos- 
sibility of fine tuning to heavy nuclei, of using EFT/RG 
to guide next-generation functional forms, and of bench- 
marking with ab-initio methods for low-momentum in- 
teractions. Work in these directions is in progress. 

In summary, we have presented the first results for nu- 
clear matter based on chiral NN and 3N interactions with 
RG evolution. The chiral EFT framework provides a sys- 
tematic improvable Hamiltonian while the softening of 
nuclear forces by RG evolution enhances the convergence 
and control of the many-body calculation. The empiri- 
cal saturation point is reproduced with theoretical uncer- 
tainties despite input only from few-body data. Because 
of the fine cancellations, however, significant reduction of 
the errors will be needed before direct DFT calculations 
of nuclei are competitive. Nevertheless, these results are 
very promising for a unified description of all nuclei and 
nuclear matter. 
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